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Trapped and cooled gases of alkali atoms can be manipulated to exhibit a variety of interesting 
phenomena. For example, dilute gases of fermionic atoms, in 2 hyperfine states, can be cooled 
to temperatures where they become superfluid. An external field can be applied to tune the scat- 
tering length a. When \a\ exceeds the interparticle spacing, nonperturbative tools are needed to 
study the system theoretically. The unitary limit, \a\ ^ oo, is particularly interesting due to its 
universality and symmetry. Lattice field theory and effective field theory can be used to system- 
atically calculate properties of this system. Results are presented for the finite temperature phase 
transition and for behavior near zero temperature. 
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1. Introduction 

In this talk 1 present results from studying a dilute gas of 2 species of nonrelativistic fermions 
interacting through a short range attraction. The most relevant experimental realization of this 
system is the trapping and cooling of fermionic isotopes of alkali atoms, a vibrant activity in many 
labs around the world. These atomic gas experiments, and the theory which describes them, can 
also be considered to be models of a dilute neutron gas. 

The Hamiltonian is the sum of kinetic energy operators and a short range potential 

^= I:V2 + £V2 +££v(|r,-r,|). (1.1) 

^ \i=l ;=l / '=17=1 

This work is concerned with equal populations of the 2 species, Ni = N2= N. We adopt nonrel- 
ativistic units h = m = I since the fermion mass m, not the speed of light, is the relevant kine- 
matic factor relating temporal and spatial units. Realistic potentials are the van der Waals potential 
v(r) ~ {ro/r)^ for atomic gases and the Yukawa potential v(r) ~ exp(— r/ro)/r for a neutron gas. 
In either case ro is the length scale characteristic of the specific potential. 

As we know, low momentum processes do not probe short distances, so under the right con- 
ditions, the many body physics is the same for atoms and neutrons. The typical momentum 
in a fermion gas is the Fermi momentum, defined through the particle density n = N /V to be 
kf = {3K^n)^^^. Quantitatively then, the details of the potential do not affect the physics when 
kpro 1, which we call the dilute regime. The atomic gases are indeed dilute; interparticle dis- 
tances are typically 2-3 orders of magnitude larger than tq. Although reaUstic neutron matter is 
never in the dilute regime - without high densities of about 0.1 GeV per (fm)^, neutrons jS -decay 
to protons - idealized dilute neutron matter is a valid theoretical limit of realistic neutron matter 
and a warm-up for nuclear matter. 

Looking at an effective range expansion of the scattering amplitude. 



— l/a + \k^rQ + ... — ik 

(defined through the S matrix as 5 = 1 + ik.s/) we can see that the only relevant parameter is 
the S-wave scattering length a in the dilute regime. (Higher partial waves are negligible at low 
energies.) Note we follow the convention is that a < corresponds to the 2-body system having no 
bound state. With — 1 <C kpa < the ground state is a Bardeen-Cooper-Schrieffer superfluid, while 
at < kpa ^ 1 tightly bound difermion molecules are the ground state and form a Bose-Einstein 
condensate at low temperatures. Mean field theory provides an accurate description in both of these 
extremes. The case of particular interest to us is ^ 1, which happens when a bound state is 
close to threshold. Here fluctuations dominate and mean field theory is of no use. In particular, 
when \a\ oo, the scattering length is no longer a physical length scale. The only length scale left 
in the problem is kp^ (or the interparticle spacing n^^^^), giving universal physics. This is called 
the unitary regime. 

The atomic physics experiments are pristine, versatile environments to study the dependence 
of many body physics on scattering length. The 2 species of fermions in the traps are 2 hyperfine 
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spin states, which react differently to an external magnetic field. By tuning this magnetic field to a 
Feshbach resonance, the scattering length can be fine-tuned to infinity. 

Large scattering lengths arise in nuclear physics as well. Low energy nucleon-nucleon scatter- 
ing is governed by one-pion exchange, so dimensional analysis would suggest that the scattering 
length should be the same size as the pion Compton wavelength, 1 .4 fm. However, in fact the nn 
scattering length is —18.5 fm, and in ?i;?-scattering = —23.76 fm and a, = 5.42 fm. Thus nature 
exhibits an order-of-magnitude of fine tuning. Direct Monte Carlo calculation from lattice QCD of 
A^A'^ scattering and binding as a function of quark mass is underway and should pin down how this 
fine tuning arises from QCD. The study of the unitary Fermi gas is relevant for nuclear physics in 
the dilute, \a\ oo, and oo limits. 

Recent experiments have cooled atomic gases in the unitary regime to temperatures near the 
superfluid transition [|], ^ ||, 0, ^. For example, striking evidence of quantized vortices signal 
superfluidity in a rotating gas of lithium-6 Here we discuss the application of techniques 
familiar to us at this conference toward this interesting system. 



2. Lattice field theory at 

Since the details of the realistic potential in ( |1.1| ) are irrelevant for dilute gases, we can re- 
place it by a local 4-fermion interaction with a tunable bare coefficient Cq. Then the discretized 
Lagrangian is the same as in the attractive 3-dimensional Hubbard model 

^ = V^^(^,-^v2-A^)^/.+ ^(VAtvA)^ (2.1) 

An attractive interaction corresponds to Co < 0. We have introduced an anisotropy relating spatial 
and temporal lattice spacings: ^ = b^/bf. An external source for pairing can be included in the 
Lagrangian: 7vAi/A+h.c. Doing so allows one to compute the condensate £ = (yy) finite volume. 
The result in the limit of V ^ oo, / ^ is an order parameter for the spontaneous breaking of 
fermion number. Although the infinite volume limit must be taken before 7 — > limit, in this work 
we assume J is not small enough to see finite volume effects. 

In order to integrate out the fermion degrees-of-freedom, an auxiliary scalar field <p is intro- 
duced to complete the square 

^ivWy^-^/-<p¥W- (2.2) 

The result is a partition function which can be expressed as a path integral with a nonnegative 
integrand: 

Z= f&(l)det[\J\^ + K^K]e-^K (2.3) 



It has always been known that the attractive Hubbard model has no sign problem, but interest lay in 
the repulsive version as a model of superconductivity in electron systems. For dilute Fermi gases, 
this theory provides a complete effective field theory with a continuum limit [^. 

The 4-fermion coupling Co can be tuned to obtain any scattering length a through 

1 t f d^p 1 



Co^ 47ia Jbz (271)3 |p|2(l + |p|2/4§) 
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= 0.14 
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= 0.155 




= 0.175 




Figure 1: The order parameter, after extrapolating 7 — > 0, as a function of anisotropy for 4 values of the 
coupling. A rapid change in E is apparent between broken (superfluid) and restored (normal) phases. 



where pj = {2/b)sm{pjb/2). The matching is obtained by requiring 2—^2 scattering calculated 
in the vacuum from ( [2.1[ ) to give the scattering amplitude ([T^. 

A year ago I performed an exploratory Monte Carlo calculation of the critical temperature 
separating the superfluid and normal phases The lattice volume and chemical potential were 
fixed to 8^ X 16 and ixb^ /E, =0.4 respectively. As a result the fermion number density was approx- 
imately 1 fermion per 4-5 spatial lattice sites. Clearly more dilute calculations are desirable to take 
the continuum limit. Nevertheless, after extrapolating the external pairing source 7 ^ 0, a clear 
jump in the pairing condensate is observed as a function of anisotropy, or temperature. Figure |] 
shows the condensate vanishing to zero across the transition as the U(l) fermion number symmetry 
is restored. The four different curves correspond to four values of the bare coupling. Since the 



matching condition (2.4) also includes the anisotropy E,, the physical scattering length is varying in 
the critical region, broadened out due to finite volume effects. In order to work at fixed scattering 
length it may be more convenient to hold both and Co fixed and vary /i across the transition. A 
back-of-the-envelope conversion of the critical region observed in Fig. [l| into a critical temperature 
yields Tc/Tp ^ 0.04 around l/\a\ < 0.5. 

In previous work, the attractive 3D Hubbard model has been studied in the condensed matter 
literature |^], and similar numerical studies have attempted to model nuclear matter | IC]. 

An advance was made recently in Monte Carlo algorithms for this system [11, 12]. Dia- 
grammatic Monte Carlo utilizes the convergence of the perturbative expansion for any value of 
the coupling, using importance sampling to generate configurations of vertices and propagators in 
position space. The worm updating algorithm in particular leads to a shortening of autocorrelation 
times [p^]. Calculations have been done using a finite scaling analysis and have taken the contin- 
uum limit. They find Tc/Tf = 0. 152(7) [|T^, |TJ. Note that at finite lattice spacings they find much 
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lower values for Tc, in rough agreement with my exploratory calculation. The extrapolation to the 
continuum limit relies on correctly fitting the sizable discretization effects. 

A calculation using different Monte Carlo methods quotes T^/Tf = 0.23(2) [|T5|], while an 
upper bound Tc/Tf < 0.14 was placed in another work [|T6|]. The systematic errors appear under 
tightest control in the work of Burovski et ai, but discrepancy between the various methods needs 
to be understood and resolved. 



3. Effective field theory at T = 



(This section presents collaborative work undertaken with D. T. Son [17].) 

Now let us consider temperatures far below the critical temperature, to where the only rel- 
evant degrees of freedom are the phonons, fluctuations (p{t,x) in the phase of the zero temper- 
ature condensate {iffY) = I (V'^'/'^)! exp(— 2/(p). Phonons dominate the physics when the thermal 
wavelength is much larger than the coherence length, roughly the size of the correlated fermion 
pair: y/ln/T » Vf/Ao, where the Fermi velocity vp is defined in terms of the fermion number 
density and Aq is the gap in the fermion spectrum. The phonon is a massless Goldstone boson 
resulting from the spontaneous breaking of the U(l) particle number symmetry: y ~^ exp(/a)i/A, 
i//^^exp(— /a). Below we construct an effective field theory for the phonons in a unitary 
Fermi gas. 

In using symmetries to constrain terms in the phonon Lagrangian, it is convenient to absorb 
the chemical potential into a field 6{t,x) = pit — (p{t,x). The lowest-order phonon Lagrangian is 

^0 = coidrd-^—^] . (3.1) 



Galilean invariance forces the first derivatives of 6 to appear in the linear combination X = dfd — 
|V6p/2. Scale invariance dictates the functional form of at the unitary point l/a = 0; dimen- 
sional analysis permits other powers of X when a is finite. The coefficient cq must be determined 
from experiment or from calculation in the microscopic theory, for example from the energy per 
fermion in the ground state. Below we write cq = 2^/^/(1571^ c^^/^), where ^ is the ratio of the 
energy per particle compared to the ideal Fermi gas energy per particle. 



The lowest-order Lagrangian ( |3.1| ) is equivalent to Landau's superfluid hydrodynamics. A 
definite power counting scheme is needed to go beyond leading order. Since 6 never appears 
undifferentiated, we start by counting all first derivatives of 6 as terms of Each further 

derivative brings in a power of momentum dp d" ~ ^{p")- 

It is convenient to gauge the U(l), introducing an external field (Ao,A). We also write the 
interaction using an auxiliary scalar field a with 2 free parameters qo and ro which can be tuned 
to obtain any scattering length, 1 /c? = in particular. The microscopic Lagrangian from which we 
start is as follows: 

1 

+ qQ\lf^^fO - -g'^diOdjO - (3.2) 
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Note that ^ is written with a nontrivial spatial metric gij. These external fields, Ao,A;,g,y, are 
useful tools for analyzing the system's symmetry, taking the limit where A,- — > and gij —>■ 5ij to 
obtain physical results. We may wish to keep Aq nonzero to include effects of the atomic trap. Aq 
always appears in combination with the chemical potential, so that /i — Ao(x) is a local effective 
chemical potential, e.g. 

X = ^-Ao{x)-dr(p-^{Di(p)\ (3.3) 

The microscopic action / dtdx^^ is invariant under nonrelativistic general coordinate and 
conformal transformations. The low energy effective action must be similarly invariant. The re- 
sulting next-to-leading order Lagrangian, after taking A,- — > and gij dij, is JItI ] 

^2 = ci + C2VX [{V^cpf - 9V^Ao] . (3.4) 

a/a 

Summarizing results obtained with the NLO superfluid Lagrangian, we find the leading cor- 
rection to a linear phonon dispersion relation 



(0 



+ 



The static density response function is found to be 



1 + 271^ 



And the static transverse response function is 



9C2 



Inq). 



^{q^lnq). 



(3.5) 



(3.6) 



(3.7) 



A dispersion relation requires that X^{q) < 0; consequently C2 must be positive. The predictivity 
of the effective field theory is apparent since the 3 observables over-constrain the 2 free parameters. 
Present experiments are not able to go to low enough temperatures to test these predictions. How- 
ever, Monte Carlo calculations should be able to bridge the gap between low temperatures, where 
phonons govern the physics, to temperatures around Tc, where experiments can be performed. Per- 
haps these calculations will be able to test the predictions of this effective theory. 

The alert reader might have noticed that the logarithmic corrections due to loop diagrams are 
suppressed compared to the NLO terms by another power of p^, in contradistinction to the chiral 
logarithms. In chiral perturbation theory, the logarithms appearing at NLO arise from interaction 
terms where the pion field appears undifferentiated, e.g. {7t • dj7[){lt • d-^TC). Such terms are not 
allowed in the (abelian) U(l) effective theory. 

Let us close by emphasizing a few features of this work. It is remarkable that imposing non- 
relativistic general coordinate invariance imposes additional constraints on the phonon Lagrangian 
beyond those dictated by Galilean invariance. For example, Galilean invariance would relax a linear 
combination of operators dictated by general coordinate invariance, so that {q) would depend 
on 3, not 2, low energy constants. This does not happen with relativistic theories: any Lorentz 
invariant term may be made general coordinate invariant quite simply. However, the nonrelativistic 
transformation is just the c ^ 00 limit of the relativistic transformation. 
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One might question the physical justification of requiring the fermion Lagrangian (^^J) to 
have this coordinate invariance - after all, the microscopic degrees-of-freedom are atoms whose 
phenomenologically successful potentials are not general coordinate invariant. However, we take 
it as a valid first principle that, since the underlying dynamics in atoms are described by QED and 
QCD, which are general coordinate invariant theories, the nonrelativistic remnant of this invariance 
persists and has predictive, refutable consequences as discussed above. 

Similarly, it is a logical possibility that the conformal invariance we impose need not be re- 
aUzed by nature. (Then (V^9)^ and V^Aq would be linearly independent in (p^.) However, the 
fact that the free Schrodinger equation is scale and conformally invariant |jT8|] hints that the unitary 
Fermi gas should also possess the same symmetry. 
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